3.10.36 \(\int \frac {\cos ^2(e+f x)}{(a+a \sin (e+f x))^{3/2} (c+d \sin (e+f x))} \, dx\) [936]

3.10.36.1 Optimal result
3.10.36.2 Mathematica [C] (verified)
3.10.36.3 Rubi [A] (verified)
3.10.36.4 Maple [A] (verified)
3.10.36.5 Fricas [B] (verification not implemented)
3.10.36.6 Sympy [F(-1)]
3.10.36.7 Maxima [F]
3.10.36.8 Giac [B] (verification not implemented)
3.10.36.9 Mupad [F(-1)]

3.10.36.1 Optimal result

Integrand size = 35, antiderivative size = 123 \[ \int \frac {\cos ^2(e+f x)}{(a+a \sin (e+f x))^{3/2} (c+d \sin (e+f x))} \, dx=-\frac {2 \sqrt {2} \text {arctanh}\left (\frac {\sqrt {a} \cos (e+f x)}{\sqrt {2} \sqrt {a+a \sin (e+f x)}}\right )}{a^{3/2} (c-d) f}+\frac {2 \sqrt {c+d} \text {arctanh}\left (\frac {\sqrt {a} \sqrt {d} \cos (e+f x)}{\sqrt {c+d} \sqrt {a+a \sin (e+f x)}}\right )}{a^{3/2} (c-d) \sqrt {d} f} \]

output
-2*arctanh(1/2*cos(f*x+e)*a^(1/2)*2^(1/2)/(a+a*sin(f*x+e))^(1/2))*2^(1/2)/ 
a^(3/2)/(c-d)/f+2*arctanh(cos(f*x+e)*a^(1/2)*d^(1/2)/(c+d)^(1/2)/(a+a*sin( 
f*x+e))^(1/2))*(c+d)^(1/2)/a^(3/2)/(c-d)/f/d^(1/2)
 
3.10.36.2 Mathematica [C] (verified)

Result contains higher order function than in optimal. Order 9 vs. order 3 in optimal.

Time = 3.08 (sec) , antiderivative size = 750, normalized size of antiderivative = 6.10 \[ \int \frac {\cos ^2(e+f x)}{(a+a \sin (e+f x))^{3/2} (c+d \sin (e+f x))} \, dx=-\frac {(-1)^{3/4} \left ((8+8 i) \sqrt {d} \text {arctanh}\left (\left (\frac {1}{2}+\frac {i}{2}\right ) (-1)^{3/4} \left (-1+\tan \left (\frac {1}{4} (e+f x)\right )\right )\right )+\sqrt [4]{-1} \text {RootSum}\left [c+4 d \text {$\#$1}+2 c \text {$\#$1}^2-4 d \text {$\#$1}^3+c \text {$\#$1}^4\&,\frac {-c \sqrt {d} \log \left (-\text {$\#$1}+\tan \left (\frac {1}{4} (e+f x)\right )\right )-d^{3/2} \log \left (-\text {$\#$1}+\tan \left (\frac {1}{4} (e+f x)\right )\right )-d \sqrt {c+d} \log \left (-\text {$\#$1}+\tan \left (\frac {1}{4} (e+f x)\right )\right )-2 c \sqrt {d} \log \left (-\text {$\#$1}+\tan \left (\frac {1}{4} (e+f x)\right )\right ) \text {$\#$1}-2 d^{3/2} \log \left (-\text {$\#$1}+\tan \left (\frac {1}{4} (e+f x)\right )\right ) \text {$\#$1}-c \sqrt {c+d} \log \left (-\text {$\#$1}+\tan \left (\frac {1}{4} (e+f x)\right )\right ) \text {$\#$1}+c \sqrt {d} \log \left (-\text {$\#$1}+\tan \left (\frac {1}{4} (e+f x)\right )\right ) \text {$\#$1}^2+d^{3/2} \log \left (-\text {$\#$1}+\tan \left (\frac {1}{4} (e+f x)\right )\right ) \text {$\#$1}^2+3 d \sqrt {c+d} \log \left (-\text {$\#$1}+\tan \left (\frac {1}{4} (e+f x)\right )\right ) \text {$\#$1}^2-c \sqrt {c+d} \log \left (-\text {$\#$1}+\tan \left (\frac {1}{4} (e+f x)\right )\right ) \text {$\#$1}^3}{-d-c \text {$\#$1}+3 d \text {$\#$1}^2-c \text {$\#$1}^3}\&\right ]+\sqrt [4]{-1} \text {RootSum}\left [c+4 d \text {$\#$1}+2 c \text {$\#$1}^2-4 d \text {$\#$1}^3+c \text {$\#$1}^4\&,\frac {-c \sqrt {d} \log \left (-\text {$\#$1}+\tan \left (\frac {1}{4} (e+f x)\right )\right )-d^{3/2} \log \left (-\text {$\#$1}+\tan \left (\frac {1}{4} (e+f x)\right )\right )+d \sqrt {c+d} \log \left (-\text {$\#$1}+\tan \left (\frac {1}{4} (e+f x)\right )\right )-2 c \sqrt {d} \log \left (-\text {$\#$1}+\tan \left (\frac {1}{4} (e+f x)\right )\right ) \text {$\#$1}-2 d^{3/2} \log \left (-\text {$\#$1}+\tan \left (\frac {1}{4} (e+f x)\right )\right ) \text {$\#$1}+c \sqrt {c+d} \log \left (-\text {$\#$1}+\tan \left (\frac {1}{4} (e+f x)\right )\right ) \text {$\#$1}+c \sqrt {d} \log \left (-\text {$\#$1}+\tan \left (\frac {1}{4} (e+f x)\right )\right ) \text {$\#$1}^2+d^{3/2} \log \left (-\text {$\#$1}+\tan \left (\frac {1}{4} (e+f x)\right )\right ) \text {$\#$1}^2-3 d \sqrt {c+d} \log \left (-\text {$\#$1}+\tan \left (\frac {1}{4} (e+f x)\right )\right ) \text {$\#$1}^2+c \sqrt {c+d} \log \left (-\text {$\#$1}+\tan \left (\frac {1}{4} (e+f x)\right )\right ) \text {$\#$1}^3}{-d-c \text {$\#$1}+3 d \text {$\#$1}^2-c \text {$\#$1}^3}\&\right ]\right ) \left (\cos \left (\frac {1}{2} (e+f x)\right )+\sin \left (\frac {1}{2} (e+f x)\right )\right )^3}{2 \sqrt {d} (-c+d) f (a (1+\sin (e+f x)))^{3/2}} \]

input
Integrate[Cos[e + f*x]^2/((a + a*Sin[e + f*x])^(3/2)*(c + d*Sin[e + f*x])) 
,x]
 
output
-1/2*((-1)^(3/4)*((8 + 8*I)*Sqrt[d]*ArcTanh[(1/2 + I/2)*(-1)^(3/4)*(-1 + T 
an[(e + f*x)/4])] + (-1)^(1/4)*RootSum[c + 4*d*#1 + 2*c*#1^2 - 4*d*#1^3 + 
c*#1^4 & , (-(c*Sqrt[d]*Log[-#1 + Tan[(e + f*x)/4]]) - d^(3/2)*Log[-#1 + T 
an[(e + f*x)/4]] - d*Sqrt[c + d]*Log[-#1 + Tan[(e + f*x)/4]] - 2*c*Sqrt[d] 
*Log[-#1 + Tan[(e + f*x)/4]]*#1 - 2*d^(3/2)*Log[-#1 + Tan[(e + f*x)/4]]*#1 
 - c*Sqrt[c + d]*Log[-#1 + Tan[(e + f*x)/4]]*#1 + c*Sqrt[d]*Log[-#1 + Tan[ 
(e + f*x)/4]]*#1^2 + d^(3/2)*Log[-#1 + Tan[(e + f*x)/4]]*#1^2 + 3*d*Sqrt[c 
 + d]*Log[-#1 + Tan[(e + f*x)/4]]*#1^2 - c*Sqrt[c + d]*Log[-#1 + Tan[(e + 
f*x)/4]]*#1^3)/(-d - c*#1 + 3*d*#1^2 - c*#1^3) & ] + (-1)^(1/4)*RootSum[c 
+ 4*d*#1 + 2*c*#1^2 - 4*d*#1^3 + c*#1^4 & , (-(c*Sqrt[d]*Log[-#1 + Tan[(e 
+ f*x)/4]]) - d^(3/2)*Log[-#1 + Tan[(e + f*x)/4]] + d*Sqrt[c + d]*Log[-#1 
+ Tan[(e + f*x)/4]] - 2*c*Sqrt[d]*Log[-#1 + Tan[(e + f*x)/4]]*#1 - 2*d^(3/ 
2)*Log[-#1 + Tan[(e + f*x)/4]]*#1 + c*Sqrt[c + d]*Log[-#1 + Tan[(e + f*x)/ 
4]]*#1 + c*Sqrt[d]*Log[-#1 + Tan[(e + f*x)/4]]*#1^2 + d^(3/2)*Log[-#1 + Ta 
n[(e + f*x)/4]]*#1^2 - 3*d*Sqrt[c + d]*Log[-#1 + Tan[(e + f*x)/4]]*#1^2 + 
c*Sqrt[c + d]*Log[-#1 + Tan[(e + f*x)/4]]*#1^3)/(-d - c*#1 + 3*d*#1^2 - c* 
#1^3) & ])*(Cos[(e + f*x)/2] + Sin[(e + f*x)/2])^3)/(Sqrt[d]*(-c + d)*f*(a 
*(1 + Sin[e + f*x]))^(3/2))
 
3.10.36.3 Rubi [A] (verified)

Time = 0.85 (sec) , antiderivative size = 127, normalized size of antiderivative = 1.03, number of steps used = 10, number of rules used = 9, \(\frac {\text {number of rules}}{\text {integrand size}}\) = 0.257, Rules used = {3042, 3395, 3042, 3464, 3042, 3128, 219, 3252, 221}

Below are the steps used by Rubi to obtain the solution. The rule number used for the transformation is given above next to the arrow. The rules definitions used are listed below.

\(\displaystyle \int \frac {\cos ^2(e+f x)}{(a \sin (e+f x)+a)^{3/2} (c+d \sin (e+f x))} \, dx\)

\(\Big \downarrow \) 3042

\(\displaystyle \int \frac {\cos (e+f x)^2}{(a \sin (e+f x)+a)^{3/2} (c+d \sin (e+f x))}dx\)

\(\Big \downarrow \) 3395

\(\displaystyle \frac {\int \frac {a-a \sin (e+f x)}{\sqrt {\sin (e+f x) a+a} (c+d \sin (e+f x))}dx}{a^2}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {\int \frac {a-a \sin (e+f x)}{\sqrt {\sin (e+f x) a+a} (c+d \sin (e+f x))}dx}{a^2}\)

\(\Big \downarrow \) 3464

\(\displaystyle \frac {\frac {2 a \int \frac {1}{\sqrt {\sin (e+f x) a+a}}dx}{c-d}-\frac {(c+d) \int \frac {\sqrt {\sin (e+f x) a+a}}{c+d \sin (e+f x)}dx}{c-d}}{a^2}\)

\(\Big \downarrow \) 3042

\(\displaystyle \frac {\frac {2 a \int \frac {1}{\sqrt {\sin (e+f x) a+a}}dx}{c-d}-\frac {(c+d) \int \frac {\sqrt {\sin (e+f x) a+a}}{c+d \sin (e+f x)}dx}{c-d}}{a^2}\)

\(\Big \downarrow \) 3128

\(\displaystyle \frac {-\frac {4 a \int \frac {1}{2 a-\frac {a^2 \cos ^2(e+f x)}{\sin (e+f x) a+a}}d\frac {a \cos (e+f x)}{\sqrt {\sin (e+f x) a+a}}}{f (c-d)}-\frac {(c+d) \int \frac {\sqrt {\sin (e+f x) a+a}}{c+d \sin (e+f x)}dx}{c-d}}{a^2}\)

\(\Big \downarrow \) 219

\(\displaystyle \frac {-\frac {(c+d) \int \frac {\sqrt {\sin (e+f x) a+a}}{c+d \sin (e+f x)}dx}{c-d}-\frac {2 \sqrt {2} \sqrt {a} \text {arctanh}\left (\frac {\sqrt {a} \cos (e+f x)}{\sqrt {2} \sqrt {a \sin (e+f x)+a}}\right )}{f (c-d)}}{a^2}\)

\(\Big \downarrow \) 3252

\(\displaystyle \frac {\frac {2 a (c+d) \int \frac {1}{a (c+d)-\frac {a^2 d \cos ^2(e+f x)}{\sin (e+f x) a+a}}d\frac {a \cos (e+f x)}{\sqrt {\sin (e+f x) a+a}}}{f (c-d)}-\frac {2 \sqrt {2} \sqrt {a} \text {arctanh}\left (\frac {\sqrt {a} \cos (e+f x)}{\sqrt {2} \sqrt {a \sin (e+f x)+a}}\right )}{f (c-d)}}{a^2}\)

\(\Big \downarrow \) 221

\(\displaystyle \frac {\frac {2 \sqrt {a} \sqrt {c+d} \text {arctanh}\left (\frac {\sqrt {a} \sqrt {d} \cos (e+f x)}{\sqrt {c+d} \sqrt {a \sin (e+f x)+a}}\right )}{\sqrt {d} f (c-d)}-\frac {2 \sqrt {2} \sqrt {a} \text {arctanh}\left (\frac {\sqrt {a} \cos (e+f x)}{\sqrt {2} \sqrt {a \sin (e+f x)+a}}\right )}{f (c-d)}}{a^2}\)

input
Int[Cos[e + f*x]^2/((a + a*Sin[e + f*x])^(3/2)*(c + d*Sin[e + f*x])),x]
 
output
((-2*Sqrt[2]*Sqrt[a]*ArcTanh[(Sqrt[a]*Cos[e + f*x])/(Sqrt[2]*Sqrt[a + a*Si 
n[e + f*x]])])/((c - d)*f) + (2*Sqrt[a]*Sqrt[c + d]*ArcTanh[(Sqrt[a]*Sqrt[ 
d]*Cos[e + f*x])/(Sqrt[c + d]*Sqrt[a + a*Sin[e + f*x]])])/((c - d)*Sqrt[d] 
*f))/a^2
 

3.10.36.3.1 Defintions of rubi rules used

rule 219
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(1/(Rt[a, 2]*Rt[-b, 2]))* 
ArcTanh[Rt[-b, 2]*(x/Rt[a, 2])], x] /; FreeQ[{a, b}, x] && NegQ[a/b] && (Gt 
Q[a, 0] || LtQ[b, 0])
 

rule 221
Int[((a_) + (b_.)*(x_)^2)^(-1), x_Symbol] :> Simp[(Rt[-a/b, 2]/a)*ArcTanh[x 
/Rt[-a/b, 2]], x] /; FreeQ[{a, b}, x] && NegQ[a/b]
 

rule 3042
Int[u_, x_Symbol] :> Int[DeactivateTrig[u, x], x] /; FunctionOfTrigOfLinear 
Q[u, x]
 

rule 3128
Int[1/Sqrt[(a_) + (b_.)*sin[(c_.) + (d_.)*(x_)]], x_Symbol] :> Simp[-2/d 
Subst[Int[1/(2*a - x^2), x], x, b*(Cos[c + d*x]/Sqrt[a + b*Sin[c + d*x]])], 
 x] /; FreeQ[{a, b, c, d}, x] && EqQ[a^2 - b^2, 0]
 

rule 3252
Int[Sqrt[(a_) + (b_.)*sin[(e_.) + (f_.)*(x_)]]/((c_.) + (d_.)*sin[(e_.) + ( 
f_.)*(x_)]), x_Symbol] :> Simp[-2*(b/f)   Subst[Int[1/(b*c + a*d - d*x^2), 
x], x, b*(Cos[e + f*x]/Sqrt[a + b*Sin[e + f*x]])], x] /; FreeQ[{a, b, c, d, 
 e, f}, x] && NeQ[b*c - a*d, 0] && EqQ[a^2 - b^2, 0] && NeQ[c^2 - d^2, 0]
 

rule 3395
Int[cos[(e_.) + (f_.)*(x_)]^2*((a_) + (b_.)*sin[(e_.) + (f_.)*(x_)])^(m_)*( 
(c_) + (d_.)*sin[(e_.) + (f_.)*(x_)])^(n_), x_Symbol] :> Simp[1/b^2   Int[( 
a + b*Sin[e + f*x])^(m + 1)*(c + d*Sin[e + f*x])^n*(a - b*Sin[e + f*x]), x] 
, x] /; FreeQ[{a, b, c, d, e, f, m, n}, x] && EqQ[a^2 - b^2, 0] && Integers 
Q[2*m, 2*n]
 

rule 3464
Int[((A_.) + (B_.)*sin[(e_.) + (f_.)*(x_)])/(Sqrt[(a_) + (b_.)*sin[(e_.) + 
(f_.)*(x_)]]*((c_.) + (d_.)*sin[(e_.) + (f_.)*(x_)])), x_Symbol] :> Simp[(A 
*b - a*B)/(b*c - a*d)   Int[1/Sqrt[a + b*Sin[e + f*x]], x], x] + Simp[(B*c 
- A*d)/(b*c - a*d)   Int[Sqrt[a + b*Sin[e + f*x]]/(c + d*Sin[e + f*x]), x], 
 x] /; FreeQ[{a, b, c, d, e, f, A, B}, x] && NeQ[b*c - a*d, 0] && EqQ[a^2 - 
 b^2, 0] && NeQ[c^2 - d^2, 0]
 
3.10.36.4 Maple [A] (verified)

Time = 0.25 (sec) , antiderivative size = 160, normalized size of antiderivative = 1.30

method result size
default \(-\frac {2 \left (1+\sin \left (f x +e \right )\right ) \sqrt {-a \left (\sin \left (f x +e \right )-1\right )}\, \left (\sqrt {2}\, \operatorname {arctanh}\left (\frac {\sqrt {-a \left (\sin \left (f x +e \right )-1\right )}\, \sqrt {2}}{2 \sqrt {a}}\right ) \sqrt {a \left (c +d \right ) d}-\operatorname {arctanh}\left (\frac {\sqrt {-a \left (\sin \left (f x +e \right )-1\right )}\, d}{\sqrt {a \left (c +d \right ) d}}\right ) \sqrt {a}\, c -\operatorname {arctanh}\left (\frac {\sqrt {-a \left (\sin \left (f x +e \right )-1\right )}\, d}{\sqrt {a \left (c +d \right ) d}}\right ) \sqrt {a}\, d \right )}{a^{\frac {3}{2}} \left (c -d \right ) \sqrt {a \left (c +d \right ) d}\, \cos \left (f x +e \right ) \sqrt {a +a \sin \left (f x +e \right )}\, f}\) \(160\)

input
int(cos(f*x+e)^2/(a+a*sin(f*x+e))^(3/2)/(c+d*sin(f*x+e)),x,method=_RETURNV 
ERBOSE)
 
output
-2/a^(3/2)*(1+sin(f*x+e))*(-a*(sin(f*x+e)-1))^(1/2)*(2^(1/2)*arctanh(1/2*( 
-a*(sin(f*x+e)-1))^(1/2)*2^(1/2)/a^(1/2))*(a*(c+d)*d)^(1/2)-arctanh((-a*(s 
in(f*x+e)-1))^(1/2)*d/(a*(c+d)*d)^(1/2))*a^(1/2)*c-arctanh((-a*(sin(f*x+e) 
-1))^(1/2)*d/(a*(c+d)*d)^(1/2))*a^(1/2)*d)/(c-d)/(a*(c+d)*d)^(1/2)/cos(f*x 
+e)/(a+a*sin(f*x+e))^(1/2)/f
 
3.10.36.5 Fricas [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 202 vs. \(2 (100) = 200\).

Time = 0.38 (sec) , antiderivative size = 667, normalized size of antiderivative = 5.42 \[ \int \frac {\cos ^2(e+f x)}{(a+a \sin (e+f x))^{3/2} (c+d \sin (e+f x))} \, dx=\left [-\frac {\sqrt {\frac {c + d}{a d}} \log \left (\frac {d^{2} \cos \left (f x + e\right )^{3} - {\left (6 \, c d + 7 \, d^{2}\right )} \cos \left (f x + e\right )^{2} - c^{2} - 2 \, c d - d^{2} - 4 \, {\left (d^{2} \cos \left (f x + e\right )^{2} - c d - 3 \, d^{2} - {\left (c d + 2 \, d^{2}\right )} \cos \left (f x + e\right ) + {\left (d^{2} \cos \left (f x + e\right ) + c d + 3 \, d^{2}\right )} \sin \left (f x + e\right )\right )} \sqrt {a \sin \left (f x + e\right ) + a} \sqrt {\frac {c + d}{a d}} - {\left (c^{2} + 8 \, c d + 9 \, d^{2}\right )} \cos \left (f x + e\right ) + {\left (d^{2} \cos \left (f x + e\right )^{2} - c^{2} - 2 \, c d - d^{2} + 2 \, {\left (3 \, c d + 4 \, d^{2}\right )} \cos \left (f x + e\right )\right )} \sin \left (f x + e\right )}{d^{2} \cos \left (f x + e\right )^{3} + {\left (2 \, c d + d^{2}\right )} \cos \left (f x + e\right )^{2} - c^{2} - 2 \, c d - d^{2} - {\left (c^{2} + d^{2}\right )} \cos \left (f x + e\right ) + {\left (d^{2} \cos \left (f x + e\right )^{2} - 2 \, c d \cos \left (f x + e\right ) - c^{2} - 2 \, c d - d^{2}\right )} \sin \left (f x + e\right )}\right ) + \frac {2 \, \sqrt {2} \log \left (-\frac {\cos \left (f x + e\right )^{2} - {\left (\cos \left (f x + e\right ) - 2\right )} \sin \left (f x + e\right ) + \frac {2 \, \sqrt {2} \sqrt {a \sin \left (f x + e\right ) + a} {\left (\cos \left (f x + e\right ) - \sin \left (f x + e\right ) + 1\right )}}{\sqrt {a}} + 3 \, \cos \left (f x + e\right ) + 2}{\cos \left (f x + e\right )^{2} - {\left (\cos \left (f x + e\right ) + 2\right )} \sin \left (f x + e\right ) - \cos \left (f x + e\right ) - 2}\right )}{\sqrt {a}}}{2 \, {\left (a c - a d\right )} f}, \frac {\sqrt {-\frac {c + d}{a d}} \arctan \left (\frac {\sqrt {a \sin \left (f x + e\right ) + a} {\left (d \sin \left (f x + e\right ) - c - 2 \, d\right )} \sqrt {-\frac {c + d}{a d}}}{2 \, {\left (c + d\right )} \cos \left (f x + e\right )}\right ) - \frac {\sqrt {2} \log \left (-\frac {\cos \left (f x + e\right )^{2} - {\left (\cos \left (f x + e\right ) - 2\right )} \sin \left (f x + e\right ) + \frac {2 \, \sqrt {2} \sqrt {a \sin \left (f x + e\right ) + a} {\left (\cos \left (f x + e\right ) - \sin \left (f x + e\right ) + 1\right )}}{\sqrt {a}} + 3 \, \cos \left (f x + e\right ) + 2}{\cos \left (f x + e\right )^{2} - {\left (\cos \left (f x + e\right ) + 2\right )} \sin \left (f x + e\right ) - \cos \left (f x + e\right ) - 2}\right )}{\sqrt {a}}}{{\left (a c - a d\right )} f}\right ] \]

input
integrate(cos(f*x+e)^2/(a+a*sin(f*x+e))^(3/2)/(c+d*sin(f*x+e)),x, algorith 
m="fricas")
 
output
[-1/2*(sqrt((c + d)/(a*d))*log((d^2*cos(f*x + e)^3 - (6*c*d + 7*d^2)*cos(f 
*x + e)^2 - c^2 - 2*c*d - d^2 - 4*(d^2*cos(f*x + e)^2 - c*d - 3*d^2 - (c*d 
 + 2*d^2)*cos(f*x + e) + (d^2*cos(f*x + e) + c*d + 3*d^2)*sin(f*x + e))*sq 
rt(a*sin(f*x + e) + a)*sqrt((c + d)/(a*d)) - (c^2 + 8*c*d + 9*d^2)*cos(f*x 
 + e) + (d^2*cos(f*x + e)^2 - c^2 - 2*c*d - d^2 + 2*(3*c*d + 4*d^2)*cos(f* 
x + e))*sin(f*x + e))/(d^2*cos(f*x + e)^3 + (2*c*d + d^2)*cos(f*x + e)^2 - 
 c^2 - 2*c*d - d^2 - (c^2 + d^2)*cos(f*x + e) + (d^2*cos(f*x + e)^2 - 2*c* 
d*cos(f*x + e) - c^2 - 2*c*d - d^2)*sin(f*x + e))) + 2*sqrt(2)*log(-(cos(f 
*x + e)^2 - (cos(f*x + e) - 2)*sin(f*x + e) + 2*sqrt(2)*sqrt(a*sin(f*x + e 
) + a)*(cos(f*x + e) - sin(f*x + e) + 1)/sqrt(a) + 3*cos(f*x + e) + 2)/(co 
s(f*x + e)^2 - (cos(f*x + e) + 2)*sin(f*x + e) - cos(f*x + e) - 2))/sqrt(a 
))/((a*c - a*d)*f), (sqrt(-(c + d)/(a*d))*arctan(1/2*sqrt(a*sin(f*x + e) + 
 a)*(d*sin(f*x + e) - c - 2*d)*sqrt(-(c + d)/(a*d))/((c + d)*cos(f*x + e)) 
) - sqrt(2)*log(-(cos(f*x + e)^2 - (cos(f*x + e) - 2)*sin(f*x + e) + 2*sqr 
t(2)*sqrt(a*sin(f*x + e) + a)*(cos(f*x + e) - sin(f*x + e) + 1)/sqrt(a) + 
3*cos(f*x + e) + 2)/(cos(f*x + e)^2 - (cos(f*x + e) + 2)*sin(f*x + e) - co 
s(f*x + e) - 2))/sqrt(a))/((a*c - a*d)*f)]
 
3.10.36.6 Sympy [F(-1)]

Timed out. \[ \int \frac {\cos ^2(e+f x)}{(a+a \sin (e+f x))^{3/2} (c+d \sin (e+f x))} \, dx=\text {Timed out} \]

input
integrate(cos(f*x+e)**2/(a+a*sin(f*x+e))**(3/2)/(c+d*sin(f*x+e)),x)
 
output
Timed out
 
3.10.36.7 Maxima [F]

\[ \int \frac {\cos ^2(e+f x)}{(a+a \sin (e+f x))^{3/2} (c+d \sin (e+f x))} \, dx=\int { \frac {\cos \left (f x + e\right )^{2}}{{\left (a \sin \left (f x + e\right ) + a\right )}^{\frac {3}{2}} {\left (d \sin \left (f x + e\right ) + c\right )}} \,d x } \]

input
integrate(cos(f*x+e)^2/(a+a*sin(f*x+e))^(3/2)/(c+d*sin(f*x+e)),x, algorith 
m="maxima")
 
output
integrate(cos(f*x + e)^2/((a*sin(f*x + e) + a)^(3/2)*(d*sin(f*x + e) + c)) 
, x)
 
3.10.36.8 Giac [B] (verification not implemented)

Leaf count of result is larger than twice the leaf count of optimal. 215 vs. \(2 (100) = 200\).

Time = 0.38 (sec) , antiderivative size = 215, normalized size of antiderivative = 1.75 \[ \int \frac {\cos ^2(e+f x)}{(a+a \sin (e+f x))^{3/2} (c+d \sin (e+f x))} \, dx=\frac {\sqrt {2} \sqrt {a} {\left (\frac {\sqrt {2} {\left (c + d\right )} \arctan \left (\frac {\sqrt {2} d \sin \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )}{\sqrt {-c d - d^{2}}}\right )}{{\left (a^{2} c \mathrm {sgn}\left (\cos \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) - a^{2} d \mathrm {sgn}\left (\cos \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right )\right )} \sqrt {-c d - d^{2}}} + \frac {\log \left (\sin \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right ) + 1\right )}{a^{2} c \mathrm {sgn}\left (\cos \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) - a^{2} d \mathrm {sgn}\left (\cos \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right )} - \frac {\log \left (-\sin \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right ) + 1\right )}{a^{2} c \mathrm {sgn}\left (\cos \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right ) - a^{2} d \mathrm {sgn}\left (\cos \left (-\frac {1}{4} \, \pi + \frac {1}{2} \, f x + \frac {1}{2} \, e\right )\right )}\right )}}{f} \]

input
integrate(cos(f*x+e)^2/(a+a*sin(f*x+e))^(3/2)/(c+d*sin(f*x+e)),x, algorith 
m="giac")
 
output
sqrt(2)*sqrt(a)*(sqrt(2)*(c + d)*arctan(sqrt(2)*d*sin(-1/4*pi + 1/2*f*x + 
1/2*e)/sqrt(-c*d - d^2))/((a^2*c*sgn(cos(-1/4*pi + 1/2*f*x + 1/2*e)) - a^2 
*d*sgn(cos(-1/4*pi + 1/2*f*x + 1/2*e)))*sqrt(-c*d - d^2)) + log(sin(-1/4*p 
i + 1/2*f*x + 1/2*e) + 1)/(a^2*c*sgn(cos(-1/4*pi + 1/2*f*x + 1/2*e)) - a^2 
*d*sgn(cos(-1/4*pi + 1/2*f*x + 1/2*e))) - log(-sin(-1/4*pi + 1/2*f*x + 1/2 
*e) + 1)/(a^2*c*sgn(cos(-1/4*pi + 1/2*f*x + 1/2*e)) - a^2*d*sgn(cos(-1/4*p 
i + 1/2*f*x + 1/2*e))))/f
 
3.10.36.9 Mupad [F(-1)]

Timed out. \[ \int \frac {\cos ^2(e+f x)}{(a+a \sin (e+f x))^{3/2} (c+d \sin (e+f x))} \, dx=\int \frac {{\cos \left (e+f\,x\right )}^2}{{\left (a+a\,\sin \left (e+f\,x\right )\right )}^{3/2}\,\left (c+d\,\sin \left (e+f\,x\right )\right )} \,d x \]

input
int(cos(e + f*x)^2/((a + a*sin(e + f*x))^(3/2)*(c + d*sin(e + f*x))),x)
 
output
int(cos(e + f*x)^2/((a + a*sin(e + f*x))^(3/2)*(c + d*sin(e + f*x))), x)